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Abstract 



We study flavor mixing and the U{1)a anomaly in QCD at zero and finite 

temperature. Using the instanton liquid model, we show that the strength of 

the anomaly is essentially unchanged near the critical temperature for chiral 

symmetry restoration. We demonstrate that nevertheless chiral symmetry 

restoration has important consequences for the ry and 77'. In particular, the 

strange and non-strange components of the r/ unmix near Tc- The anomaly 

does not affect the strange eta, so we expect a light purely strange pseu- 

doscalar near the phase transition. 
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1. In connection with the ongoing heavy ion program at AGS and CERN it is of great 
interest to identify possible changes in hadronic properties as matter is heated up and reaches 
the critical temperature for chiral symmetry restoration. Recently, a number of authors have 
argued that such changes might be very dramatic in the rj — rj' sector (see for 
earlier work on the subject). At zero temperature, the 77' — tt mass splitting, which is (mostly) 
due to the anomaly, is larger than any other mesonic mass splitting. This means that any 
tendency towards (partial) U{1)a restoration might lead to physical effects that are more 
easily observed than changes in, for example, the p meson or nucleon channels. 

In this letter, we wish to study the axial anomaly and the rj — rj' system at finite tem- 
perature in the instanton liquid model [p|-[TT[|. Since instantons provide the mechanism for 
the U{l)ji anomaly in QCD, the applicability of the model appears obvious. Nevertheless, 
we would like to make two additional comments. First, the model not only accounts for 
the anomaly, it gives a very successful description of hadronic phenomenology in general. 
In addition to that, the model describes spontaneous chiral symmetry breaking as well as 
its restoration at a critical temperature ~ 140 MeV. Our second comment concerns the 
anomaly. A quantitative description of the rj — rj' system in the instanton model is more 
subtle than one might expect. Indeed, the simplest, random, instanton model fails to give 
an acceptable description of the 77' mass and t] — rj' mixing. 

The plan of this letter is as follows. First, we discuss the rj—r]' system at zero temperature. 
Next, we review general arguments concerning the anomaly at finite temperature. Finally, 
we study 77 and rj' correlation functions at finite temperature. 

2. Before we go into detail, we should remind the reader of 't Hooft's mechanism for 
U{1)a breaking in QCD [|1^. The QCD partition function receives contributions from special 



field configurations, instantons, that carry topological charge. In the field of an instanton, 
the Dirac operator has a chiral zero mode, ip4'o = 0. In the case of an instanton, the 
zero mode is left handed 750o = ~4'o, while in the case of an anti-instanton it is right 
handed. The presence of a zero mode implies that for massless fermions, the amplitude 
for an isolated instanton vanishes, because it contains a factor det{ip). However, when we 



calculate a U{1)a violating observable, like the expectation value of the 't Hooft operator 
Odet = detf^tpLipR) + {L R), the Nf zero modes in the determinant cancel against Nf 
zero modes in the quark propagators. As a result, the instanton can absorb A'^^ left handed 
quarks and turn them into right handed quarks, violating axial charge by 2Nf units. 

In order to study the effect of the anomaly on the spectrum of pseudoscalar mesons, we 
have to consider correlation functions of the SU{3) singlet and octet pseudoscalar meson 
currents jo,8 = ^TsAo.s?- The diagonal singlet and octet, as well as the off-diagonal singlet- 
octet mixing correlators are given by 

Hoo = J {2ni^ + n^„„ - Aui, - n^,, - 4n- } , (i) 
iiss= I {2n:f,„ + 4n:„„ - 4nL - 4n^.. - sn- } , (2) 
nos = ^ {2n«„„ - 2ui^ - Aui, + 2n^,, + 2n:^^j , (3) 

where n4,„ = (Tr {x,y)'^^Sf {y,x)'^^ ) is a connected correlation function and Il^f^ 



Tr 



^■^(x, x)75 Ti[S^{y,y)'^^j is a disconnected correlator. Here, S^{x,y) is the quark 
propagator of a quark with flavor / and (.) denotes averaging over all gauge field configura- 
tions. In deriving we have assumed exact isospin symmetry. For comparison, the pion 
correlator is given by IItt = 11"^^ and the kaon correlator by Hk = I^con- Also, the correlators 
of the strange and non-strange components of the rj are given by IIjv^ = 211^^,^ — 11^.^ and 
= Hcon ~ ^disi while their mixing is determined by Hat^^s = A/2n^j?,. At zero temper- 
ature, we will exclusively focus on euclidean space correlation functions. This means that 
the long time behavior of the correlators is given by 11 ~ exp(— mx), where m is the ground 
state mass in the given channel. 

Due to flavor SU (3) symmetry breaking, the singlet and octet correlation functions mix. 
The physical rj and t]' states couple to linear combinations of the singlet and octet currents, 

rj' = cos 9riQ + sin Orjs. (4) 

Experimentally, the value of the mixing angle is ^ ~ —(10 — 20)°. The sign of the mixing 
angle corresponds to a reduction of the strange component of the t] and a strangeness 



enhancement in the rj . The uncertainty in the mixing angle is hard to assess. In fact, the 
concept of a mixing angle may not be very well defined, since the two states are rather far 
apart in mass. 

Qualitatively, instanton effects can be understood from the propagator in the field of a 
single instanton, S-^{x,y) = (pQ{x)(p\{y) /rrif + SNZM{.x,y) + Sm{x,y). Here, 4>o{x) is the zero 
mode wave function, Snzm{x, y) is the non-zero mode propagator and ^^(a;, y) includes mass 
corrections |T^. Let us first consider the case of exact SU{?)) symmetry. It is important 



to note that zero modes give identical contributions to the connected and disconnected 
correlators. We find IIoo = — 2 ■ /, IIos = / and IIos = 0, where we have defined 

dp<Pl f d'z\ct>{x-z)\'\ct>{y-z)\\ (5) 



Here, n{p) is the single instanton density |I4[ and z is the instanton position. Note that I is 
not singular in the chiral limit m — 0, because the instanton density contains a factor iv?. 
Strictly speaking, n{p) is proportional to m'^, but if chiral symmetry is broken, zero modes 
can be absorbed by the quark condensate. This means that we can replace one power of m 
by the effective mass m* = |'/rp^(gg) which is finite in the chiral limit. 

If we include SU{3) fiavor breaking, the situation becomes more complicated. In effective 
models, it is usually assumed that the 't Hooft interaction is SU{3) symmetric and the 
mixing is caused by quark or meson mass terms. In that case it is clear that the mixing 
angle is negative, since the mass term drives the system toward ideal mixing, corresponding 
to = —54.7°. In general, the situation is more complicated and SU (3) fiavor breaking in the 
't Hooft interaction is substantial. In the single instanton approximation, fiavor symmetry 
breaking can be taken into account by replacing the effective mass m* by m* + mf ||15[. The 



off-diagonal correlators are completely determined by disconnected contributions, since the 
zero mode contributions to the connected and disconnected correlators with the same flavor 



cancel each other. We flnd 15 



„ \/2m,, — m, ^ „ ^m* + 171^1 ^ , , 

no8 = ^ — l—^i, ^Ns,s = V2 " r, 6 



where we have replaced the current quark masses by effective masses in I*. In the single 
instanton approximation IIos is negative, which corresponds to a positive mixing angleQ. 
On the other hand, mass corrections in the non-zero mode part of the propagator give a 
positive contribution to IIos. This means that the sign of the mixing angle is determined 
by the competition between the flavor breaking in direct mass insertions and the 't Hooft 
interaction. 

In the following, we study this problem in the instanton liquid model. We refer the 
reader to ||TT|Jl6[| for details of the model. For our purposes here it is only important that 



there are three different instanton ensembles, the random (RILM), quenched (QILM) and 
unquenched (IILM) models. In the random and quenched ensembles, the topological suscep- 
tibility is finite Xtop — (^/^)) where (N/V) is the density of instantons. In the unquenched 
ensemble, topological charge is screened and Xtop ~ m{qq). Eta meson correlation functions 
in the random model[| and two different unquenched ensembles (the "stream line" and "ra- 
tio ansatz" ensembles, see |jlT|]) are shown in figure]^. We note that the U{1)a anomaly is 
"over-explained" in the random ensemble. The rj' channel is so repulsive that the correlation 
function becomes unphysical (n(x) < 0). Also, flavor symmetry breaking is very strong and 
the singlet-octet correlator is large and negative (corresponding to a positive mixing angle). 

This situation is improved in the unquenched ensembles. The fiavor-singlet correlation 
function is physical and singlet-octet mixing is smaller. The prediction for the masses and 
mixing angles depends sensitively on details of the interaction. In the streamline ensemble, 
the 7]' is still too heavy, m^> ~ 2 GeV. The mass of the rj is given by = (0.66 ± 0.12) GeV 



^ The off-diagonal correlator in the single instanton approximation was first calculated in |15|, 
but the conclusions concerning the mixing angle appear to be wrong. 

^The results in the quenched ensemble are very similar. Eta meson correlation functions in the 



random ensemble were first studied in |17|. Unfortunately, this work contains an error in the flavor 



octet and the off-diagonal singlet-octet correlation functions. In particular, IIos has the wrong sign. 
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and the mixing angle is small, ^ = (1 ± 3)°. A small mixing angle was also found in the 
Nambu and Jona-Lasinio model |T^ (see also p!9| , |20| ). In particular, these authors showed 
that a small mixing angle is not necessarily incompatible with the observed r/ — > 27 rate. 

3. A global measure of the strength of the anomaly is provided by the expectation value 
of the 't Hooft operator. In figure ^ we show the temperature dependence of {Odet) in the 
finite temperature instanton ensembles^ obtained in [|ll|] . For comparison, we also show the 
quark condensate and two different four fermion operators. At zero temperature both the 
SU{2) X SU{2) chiral and the axial U{1)a symmetry are broken. Chiral symmetry breaking 
is caused by interactions between the zero modes associated with individual instantons. As 
a result, some of the lowest states become collective and form a condensate. Near T = 125 
MeV, chiral symmetry is restored and the quark condensate goes to zero. This transition 
is due to a rearrangement of the instanton liquid, going from a disordered, random, system 
to an ensemble of topologically neutral instanton-anti-instanton pairs. Clearly, there is no 
tendency towards U{1)a restoration as chiral symmetry is restored. In fact, {Odet) has a 
maximum near the phase transition (although the uncertainty is also largest near T^). The 
reason why {Odet) survives the chiral phase transition should be clear from our discussion 
above: the 't Hooft operator can induce a tunneling event all by itsellQ. At temperatures 
significantly above the semi-classical tunneling amplitude contains the suppression factor 



n(p) ~ exp(— (2A''c/3 + A^//3)(7rpT) ) |j23| and {Odet) becomes small. This suppression factor 



is mostly due to Debye screening of the instanton field. Therefore, it does not affect the 



^ All of these ensembles have total topological charge Q = Q. This means that the topological 
susceptibility evaluated for the entire volume is not correct. This should not affect local observables. 
Indeed, we have checked that topological charge fluctuations in a sub-volume have the expected 



dependence on the quark mass and volume |21|. 
^This can also be verified in the mean field approximation to the instanton liquid for two flavors, 



see 
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instanton density below the phase transition. This was checked exphcitly by performing a 
soft pion calculation of the instanton density at small temperature and in (quenched) 
lattice calculations of the instanton density at finite temperature p5|-p7|. 



Local four fermion operators like Odet are hard to measure on the lattice. There- 
fore, a number of authors have studied f/A(l) violating mesonic susceptibilities. Here, 
mesonic susceptibilities are defined as integrals of the corresponding correlation function, 
Xr = J d'^xUrix). The most natural candidate for a U{1)a order parameter is the difference 
Xtt — Xv'y but this quantity also involves disconnected quark loops. A better observable is 
Xtt — Xs, originally suggested by the Columbia group p8|. For Nf = 2 and if chiral symmetry 
is restored, this quantity is a measure of U{1)a breaking. A nice feature of Xn — Xs is that 
it can be expressed in terms of the spectral density p(A) of the Dirac operator 



For comparison, the quark condensate is given by 



{qq) = -2m / ^A-^^^. (8) 

These results allow us to constrain the low virtuality part of the spectrum. For U{1)a to be 
broken but chiral symmetry restored, we require Xtt~Xs to be finite in the limit m — while 
(qq) goes to zero. This requirement is clearly satisfied by p(A) ~ m^f6{X), corresponding to 
a dilute system of instantons. However, interactions among zero and non-zero modes might 
alter the shape of the spectrum. It is therefore interesting to note that the criterion given 
above is also satisfied by a non-analytic spectral density p(A) ~ A" with a < 1. 

In order to study this question in more detail, we have determined the spectrum of 
the Dirac operator in the instanton liquid (for Nf = 2) for several different values of the 
temperature and the quark masses, see figure ^. Above we clearly observe a peak in 
the spectrum near A = 0. The number of eigenvalues in the peak is nicely consistent with 
A^(A ^ 0) ~ m^. Below Tc, most of the small eigenvalues are related to chiral symmetry 
breaking. Their number is proportional to the effective mass m* and almost independent of 



the current mass m. There is a dip in the spectrum for small quark masses. This is a finite 
volume effect. In a finite volume, the spectral density near zero will always go to zero as 
m —>■ 0. We have checked that the width of the dip in the spectrum decreases as the volume 
is increased. 

A number of groups have measured Xn ~ Xs (or the corresponding screening masses) 
on the lattice [p^-|30i. Most of the published results indicate that U{1)a remains broken, 



although recent results by the Columbia group have questioned that conclusion [^. From 
the literature, it is not clear whether lattice simulations find the peak in the spectrum 
observed in the instanton liquid. The Columbia group has measured the valence mass 



dependence of the quark condensate, which is a folded version of the Dirac spectrum p8 
Their result looks very smooth, not indicative of a small virtuality peak. One should note, 
however, that instanton calculations focus exclusively on the small virtuality part of the 
spectrum while the number of eigenvalues in lattice simulations is much larger. Also, both 
instanton and lattice simulations suffer from certain artefacts if the quark mass is made too 
small. In our case, if the mass is too small, isolated instantons are rare and the constraint 
Q = affects the results. On the lattice, for small quark masses one might run into problems 
with chiral fermions. 

4. The U{1)a anomaly at finite temperature is usually discussed in terms of the effective 
lagrangian |^ 

C = ^Tr ((9^<I>)(a/,$t)) - Tr + <l>t)) + V{<l><l>^) + c (det $ + det <l>^) , (9) 

where $ is a meson field in the (3, 3) representation of U{3) x f/(3), V^($<l>'^) is a U{3) x U{3) 
symmetric potential (usually taken be quartic), Ai is a mass matrix and c controls the 
strength of the U{1)a breaking interaction. If the coupling is taken to be c = Xtop/ (12/^), the 
effective lagrangian reproduces the Witten-Veneziano relation f^m^, = Xtop- In a quenched 
ensemble, we can further identify Xtop — (N/V). The temperature dependence of c is 
usually estimated from the semi-classical tunneling amplitude n(p) ~ exp(— (8/3)(7rpT)^). 
As a result, the strength of the anomaly is reduced by a factor ~ 5 near Tc. If the anomaly 
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becomes weaker, the eigenstates are determined by the mass matrix. In that case, the mixing 
angle is close to ideal 9 = —54.7° and the non-strange 1] is almost degenerate with the pion. 

There are several points in this line of argument that are not entirely correct. The 
strength of the 't Hooft term is not controlled by the topological susceptibility {xtop = in 
full QCD!), 

Xtop is not proportional to the instanton density (for the same reason), and, at 
least below Tc, the semi-classical estimate for the instanton density is not applicable. As we 
saw above, at Tc the instanton liquid is rearranged but the strength of the U{1)a anomaly 
does not change very much. However, chiral symmetry restoration affects the structure of 
flavor mixing in the 77 — 77' system (see figure The mixing between the strange and non- 
strange eta is controlled by the light quark condensate (see equation (D), so ?77vs and rjs do 
not mix above T^. As a result, the mixing angle is not close to zero, as it is at T = 0, but 
close to ideal. Furthermore, the anomaly can only affect the non-strange rj, not the strange 
one. Therefore, if the anomaly is sufficiently strong, the rj^s will be heavier than the 775. 

We should compare this scenario to other possibilities discussed in the recent literature. 
A number of authors have noticed that for three massless flavors, both the rj and rj' are 
unaffected by the anomaly above JT^p2| -p^. This is a related effect, but not that relevant 



for QCD, where the strange mass is not small. Also, two recent papers have argued that 
the ?7Ar5' and rjs decouple near Tc However, in these works the effect is caused by the 

disappearance of the anomaly and as a result, the rjs is always predicted to be heavier than 
the rjiys- The scenario proposed here is consistent the effective lagrangianQ (H). However, 
most authors employ a "first order" treatment of flavor symmetry breaking and neglect 
terms of order (m^ — m„)c. These terms are precisely what causes the effect discussed here. 
5. To explore this phenomenon in a more quantitative way, we study t] — rj' correla- 



^ Our scenario cannot be described in terms of the non-linear effective lagrangian employed in 
This should not be suprising; non-linear effective meson theories have to be used with care near 
the chiral phase transition. 
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tion functions in the instanton model at finite temperature. Tlie results were obtained for 
three flavors with masses m„ = = 22 MeV and = 155 MeV. We consider temporal 
correlation functions, rather than the spacelike screening correlators usually calculated on 
the lattice. Temporal correlators have the advantage that they are directly related to the 
spectrum of physical excitations. 

The results are shown in figure ^. Correlation functions below Tc are shown by the open 
squares (T = 55 MeV), pentagons (T = 78 MeV) and hexagons (T = 111 MeV), while the 
correlators near and above Tc are denoted by closed squares (T = 126 MeV) and pentagons 
(T = 145 MeV). Below Tc the singlet correlation function is strongly repulsive, while the 
octet correlator shows some attraction at larger distance. The off-diagonal correlator is small 
and positive, corresponding to a negative mixing angle. The strange and non-strange eta 
correlation functions are very similar, which is a sign for strong flavor mixing. This is also 
seen directly from the off-diagonal correlator between r]s and rjNs- 

Above Tc, the picture changes. The off-diagonal singlet-octet changes sign and its value 
at intermediate distances r ~ 0.5 fm is significantly larger. The strange and non-strange 
eta correlators are very different from each other. The non-strange correlation function 
is very repulsive, while the strange one is significantly larger. This clearly supports the 
scenario presented above. Near Tc the eigenstates are essentially the strange and non- 
strange components of the i], with the rjs being the lighter of the two states. This picture 
is not realized completely, IIs^Ars does not vanish and the singlet eta is still somewhat more 
repulsive than the octet eta correlation function. This is due to the fact that the light quark 
mass does not vanish. In particular, in this simulation the ratio (m^ + m^)/(2ms) = 1/7, 
which is about three times larger than the physical mass ratio. 

It is difficult to provide a quantitative analysis of temporal correlation functions in the 
vicinity of the phase transition. At high temperature the temporal direction in a euclidean 
box becomes short and there is no unique way to separate out the contribution from ex- 
cited states. Nevertheless, under some simplifying assumptions one can try to translate 
the correlation functions shown in figure |^ into definite predictions concerning the masses 
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of the T) and rj'. For definiteness, we will use ideal mixing above Tc and fix the threshold 
for the perturbative continuum at 1 GeV. In this case, the masses of the strange and non- 
strange components of the r] at T = 126 MeV are given by m^^ = (0.420 ± 0.120) GeV and 
mr,s = 1.250 ± 0.400) GeV. 

6. In summary, we studied flavor mixing and the axial anomaly at zero and finite 
temperature. At zero temperature, we have emphasized two points. The first one is that the 
correlations among instantons that lead to topological charge screening are also important 
in reproducing the rj and rj' correlation functions. The second one is that flavor symmetry 
breaking in the instanton induced interaction is substantial. It acts against the symmetry 
breaking from direct mass insertions. As a result, the rj — rj' mixing angle is small. 

At finite temperature, we have shown that the strength of the anomaly as measured by 
the expectation value of the 't Hooft operator is essentially independent of temperature below 
Tc. In our model, (Odet) even peaks near Tc and then drops at larger temperatures. In terms 
of the spectral density of the dirac operator, the anomaly is due to a spike p(A) ~ m^fS{X) 
in the spectrum at zero virtuality. 

Our most important result concerns the role of flavor mixing at finite temperature. 
Although the strength of the anomaly is not reduced near T^ chiral symmetry restoration 
affects the rj—rj' system. If the light quark condensate vanishes, transitions between light and 
strange pseudoscalars are suppressed. As a result, the eigenstates are given by the strange 
and non-strange components of the rj. The anomaly does not affect the strange r^, so we 
predict a purely strange, light pseudoscalar near the transition. We have estimated the mass 
of this state to be around 400 MeV. A light strange meson is of interest in connection with 
strangeness production in relativistic heavy ion collisions. This suggests that the coupling 
of a light r]s to kaons and etas should be studied in more detail. 

Acknowledgements: I would like to thank V. Koch and E. Shuryak for useful discussions. 
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FIGURES 

FIG. 1. Singlet and octet eta correlation functions in different instanton ensembles at zero 
temperature. All correlation functions are normalized to free quark propagation. 

FIG. 2. Temperature dependence of the quark condensate, the 't Hooft operator and two 
different four fermion operators in the instanton liquid. All condensates are normalized to their 
T = values. 

FIG. 3. Dirac spectra below and above the chiral phase transition for different dynamical quark 
masses. The spectral density is given in arbitrary units. The instanton density was held fixed at 

{N/V)A~'^ = 1, where A is the QCD scale parameter. Quark masses and inverse temperature 
P = are given in units of A. 

FIG. 4. Leading contributions to flavor mixing in the 'q — 'q' system below and above the chiral 
phase transition. 

FIG. 5. Correlation functions for singlet and octet, off-diagonal singlet-octet as well as strange 
and non-strange etas at different temperatures. Open squares, pentagons and hexagons are below 
Tc, closed squares and pentagons above T^. 
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